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ABSTRACT 

We use the publicly available subhalo catalogs from the via Lactea II simulation along with a Gpc-scale 
N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population 
of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us 
to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we 
identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified 
mass threshold before it is reionized; after this time the photoionizing background will suppress both the 
formation of stars and the accretion of gas. We show that varying the reionization time over the range expected 
for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. 
This conclusion is in contradiction with a number of studi es in t he literature, and we conclude that this is 
a result of inconsistent application of the results of iGnedinI (120001) : subtle changes in the assumptions about 
how reionization affects star formation in small galaxies can lead to large changes in the effect of changing 
the reionization time on the number of satellites. We compare our satellite galaxies to observations using 
both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and 
account for observational completeness effects. Additionally, if we assume that the mass threshold is set 
by the virial temperature Ty^y- = 8 x lO-'K we find that our model accurately matches the v^ax distribution, 
radial distribution, and luminosity function of observed Milky Way satellites for a reionization time Zreion = 8^2j 
assuming that the via Lactea II subhalo distribution is representative of the Milky Way. This results in the 
presence of 540;!;34q satellite galaxies. 

Subject headings: cosmology:theory — large-scale structure of universe — dark matter 



1. INTRODUCTION 

While the cold dark matter (CDM) paradigm has been very 
successful in explaining the large scale distribution of matter 
in the universe, one final test lies in its ability to predict the 
distribution of matter on small scales, including the distribu- 
tion of satellite galaxies around the Milky Way. In the hier- 
archical model of structure formation, massive objects such 
as the Milky Way halo are built up through a series of merg- 
ers where small, dense objects collapse early and merge to 
form larger objects. High-resolution numerical simulations 
have shown that the dense cores from a significant number 
of these small building blocks should survive today as grav- 
itationally bound subhalos (Klvpin et al.lll999 iMoore et al.l 



significant e ffect on the ability o f a subhalo to host a lumi- 
nous galaxv (iBullock et alJl2000t iGn edin 20 00t iBenson et al] 



2002'; 'Somerville '2002; Dekel & Woo 2003; Shaviv & DekeJl 
2003; Li et al. 2008; Madau et al. 2008; Koposov et al.^200i 



[1999; Diemand et al.ii2008: Springel et al.ll2008h. F rom these 
initial results jKlvpin et al.ll999l;lMoore et al.ll999h and from 
semi-analvtic modeling jlCauffmann et al.lll993b it was also 
clear that there are far fewer luminous dwarf satellites around 
the Milky Way than bound dark matter subhalos in the simu- 
lations. 

There has been no shortage of solutions proposed to ex- 
plain this discrepancy, including those which modify the 
dark matter power spectrum to reduce the amount of small 
scale power through warm dark matter-like models (e.g., 
Colin et al. 2000; Bode et al. 2001; Zentner & Bullock 2003; 
Busha et al.ll2007l) . Less exotic models, however, such as the 
presence of a number of "dark-dark halos," dark matter sub- 
halos that do not host galaxies, provide a well-motivated res- 
olution to this issue within the CDM paradigm. In particular, 
the presence of photoionizing radiation is expected to have a 



Star formation in dark matter halos in the early universe is 
only expected to be able to begin once the halo has grown 
massive enough to cool efficiently by atomic cooling, typ- 
ically around Tyir ~ lO'* K. However, before many halos 
are able to reach this mass, the universe enters the phase of 
reionization, in which photoionizing UV radiation is released 
by the early generations of stars and quasars. This ioniz- 
ing radiation heats the halo gas to a temperature of a few 
times lO'^K, preventing it from being pulled into the shallow 
potential wells of halos with virial temperatures lower than 
- 10^ K (iThoul & Weinberg|[T996t Ibiikstra et aUlIOOl , ef- 
fectively suppressing further star formation. 

Recent analysis of the Sloan Digital Sky Survey (SDSS) 
has resulted in the discovery of a large number of low- 
surface brig htness dwarf galaxies ( Willman et al. 20 05 aig 
Belokurov et al. 2006; Zucker et al. 2006a b; Belokuro v et al. | 
llrwinet al, .2007; .Walsh et al.. ,2007; .Belokurov et alj 
which is now beginning to shed a new light on the 
mapping of galaxies onto dark matter halos at the low- 
luminosity end. Several of these newly-discovered satel- 
lites have luminosities similar to those of the least lumi- 
nous globular clusters, and a dynamical analysis indicates that 
they have the largest mass-to-light ratio of any known galax- 




ies ([Ma rtin et al. 
[lOOa lGe ha e t ak 



ISimon & Gehal 120071; IStrigariet al.) 
Given both the magnitude limit and 
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sky-coverage fraction of the SDSS survey, it is certainly rea- 
sonable to assume that we have only detected a fraction of 
the M i Ucy Way s atellites ( iKoposov et al.ll2008l: iTollerud et"!!] 
I2008t IWalsh et a l. 2009), and there are exiting prospects 
for disco very of more satellites in future dee p and wide 
surve ys jAbbott et al.l|2005; Keller et a l. 2007; Jvezic et alj 
120081) . A full understanding of the mapping between lumi- 
nous satellite galaxies and dark matter subhalos will require 
a meas urement of the luminosity distribution, rad ial distri- 
bution jKravtsov et al.ll2004bt [Willman et al."2004'). and the 
kinematic properties of the satellites (Strigari et al. 200^. 
Recently, it has become understood that the universe 



Sokasian et al 
2006t iLidz et al. 





2003 




2007 



large scales (10 Mpc-1 Gpc; e . g., [S 

Barkana & Loeb|_|2004|; llUevet al.1 1200 

Alvarez et al. 200811 Given constraints on the global reioniza 
tion history, e. g. from the cosmic m icrowave background po- 
larization (e.g.T lSpergel et al.]l200€l) . there is still a significant 
uncertainty in the precise reionization epoch o f the matter in 
the Milky Way. In particular, the calculations of I Alvarez et alj 
(|2008i) indicate that there is substantial scatter in the reioniza- 
tion histories of halos of a given mass, and that on average 
Milky Way mass halos reionization redshifts have approxi- 
mately the same distribution as that of the universe as a whole. 
A change in the reionization time of the Milky Way may have 
a dramatic impact on its satellite population because the reion- 
ization history may directly affect the ability of a subhalo to 
reach the size where gas is able to cool and begin the star 
formation process. The primary aim of this paper is to un- 
derstand what, if any, effect the reionization epoch of a given 
Milky Way halo has on its satellite population. 

The effect of the reionization epoch on the satellite pop- 
ulation has been previously addressed in the literatur e , with 
most studies finding little effect (e.g., iSomervilT^ I2002t 
l&ravtsov et al. 2004t^). These studies primarily have ad- 
dressed the luminosity range of the "classical dwarfs" in the 
Milky Way. Given the dramatically different observational 
picture that has emerged with new observations from SDSS, 
combined with the possible importance of a spread in reion- 
ization epochs expected from inhomogeneous reionization, 
we re-investigate this question here. We combine a high- 
resolu tion dark matter simulation (the via Lactea II simula- 
tion of lDiemand et al.l2007bt l2008l) with various assumptions 
about star formation in small halos, and compare to up-to-date 
constraints from the full observed satellite galaxy population. 
In this work, we critically examine the assumptions about the 
rate at which photoionizing UV radiation is able to heat halo 
gas. We find that the exact time a halo reionizes can have a 
significant impact on the satellite population, and use compar- 
isons with the Milky Way's satellite distribution to constrain 
the reionization time of our own halo. 

In section §2 we discuss our simulation and models, includ- 
ing methods for identifying subhalos that host satellite galax- 
ies, and determining magnitudes of these galaxies. In §3 we 
discuss the observational sample of satellite galaxies that we 
compare to our model, including measurements and correc- 
tions for the luminosity, Vmax and radial distribution functions. 
In §4 we directly compare our modeled and observed samples 
in a manner that fairly accounts for the incompleteness of the 
observations. In §5 we discuss how our results compare with 
previous work in the literature, in particular addressing the 
differences between our results and those that have found lit- 
tle change in the satellite population for varying reionization 
times (,Somerville..2002; .Kravtsov et al...2004bi) . Finally, we 



present our conclusions in §6. 

2. SIMULATIONS AND MODELING 

2.1. Modeling the Large-scale Structure of Reionization 

In order to understand the distribution of reionization times 
of Milky Way sized galaxies, including any effect this may 
have on the satellite galaxy population, it is necessary to un- 
derstand the distribution o f dark matter on the largest scales. 
We use the recent results of lAlvarez et alj ( 120081) . in which the 
reionization process was modeled using an N-body simulation 
of a 1 Gpc box combined with an analytic prescription for 
predicting the reionization time for all points in the box, as de- 
scribed below. For more det ails on the reionizatio n simulation 
and the halo correlation, see lAlvarez et all (l2008b . 

The N-body simulation used the code Gadget2 jSpringell 
I2005h to evolve 1120'' dark matter particles in a cosmol- 
ogy with fi„, = 0.25, JIa = 0.75, and = 0.8, with particle 
mass resolution Mp = 4.94 x 1O'"/!"'M0. The initial condi- 
tions were generat ed using the 2nd or der Lagrangian pertur- 
bation code 2LPT dCrocce et al.ll2006h . A Frie nds-of-Friends 
group finder based on the Ntropy framework ([Gardner et alj 
120071) was run on the z = output and identified all halos with 
MfoF > 1 .58 X lO'^/i"'M0 = 32 particles. This simulation was 
run in conjunction with the LasDamas collaboration'. The 
reionization history for this simulation was then calculated us- 
ing the density field of the initial conditions. Working at a sin- 
gle point in space, the dark matter density field is smoothed 
over a series increasing radii and use the EPS formalism to 
ask, for a given redshift, what is the smallest (if any) radius 
at which the fraction of collapsed mass in halos greater than 
some threshold is greater than a specified fraction, 1 /C-. The 
reionization time, Zreion, for that point is the earliest redshift at 
which this criteria is first met for any radius. For this simula- 
tion we use C,c= 10 and a threshold mass IO^Mq. Combined 
with the group catalog, this results in a reionization time for 
each halo in the simulation. 

2.2. Modeling Galaxy Formation 

While the above simulations are sufficient for measuring 
the reionization history of Milky Way mass halos, signif- 
icantly higher mass resolution is needed to understand the 
subhalo distribution. For this, we use the publicly avail- 
able mass accretion h istories^ from the hig h- resolut ion via 
Lactea II simulation dPiemand et alj 12007 ailbl l2008h . This 
simulation models a single dark matter halo with virial mass 
Miudo = 1 .9 X lO'^/r'M© that is able to resolve subhalos down 
to mass limit M,,,/, > 10''/!~'Mq. The publicly available data 
includes the distribution, tidal mass, and Vmax histories for the 
most massive progenitors of all z = subhalos back to redshift 
z = 28. 

In order to connect the via Lactea II dark matter sub- 
halo population to a satellite galaxy population, we assume 
that stars begin forming once atomic cooling becomes ef- 
fective, when a halo shock heats to virial temperature k, 
8 X \(fiK, but that reionization heats the g as in the subhalos to 
the point where this becomes ineffective (|Thoul & Weinberg! 
[1991 iKepneretanfTggI IWise & Abell[2008^ . We treat this 
heating as an instantaneous process (see section 5 for a dis- 
cussion), causing reionization to end star formation for the 

' http://lss.phy.vanderbilt.edu/lasdamas/ 

^ http://www.ucolick.org/'^diemand/vl/data.html 
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vast majority of the ^ 2500 potential satellite galaxies; for 
these low-mass halos, all star formation must happen before 
Zieion- With this in mind, we can define a subhalo as being a 
satellite galaxy using a two parameter model: A subhalo must 
grow to a threshold mass, Mt, above which HI cooling will 
allow star formation, before the host halo reionizes at Zreion in 
order to host a satellite. 

While we demonstrate the effects of varying both p aram- 
eters in the next section, the work of I Abel et al.l (120021) uses 
high resolution AMR simulations to model the formation of 
the first stars and indicates that we anticipate Mt « 10^ - 
10^/i~'Mo. It is important to note that this process of hy- 
drogen cooling simply defines a minimum mass of the pop- 
ulation of the dark matter subhalos that could host satellite 
galaxies. However, this work predicts the stars forming in 
these halos to be very massive and short-lived. As such 
these very first star forming halos cannot be the direct pro- 
genitors of Milky Way satellites, which are observed to be 
metal-enriched objects with stars presumably of masses less 
than a solar mass. More relevant here are the calculations of 
|Wise & Abel (2008), who followed the build up of halos up 
to the masses when they start cooling via Lyman-alpha from 
neutral hydrogen. They included the radiative as well as the 
supernova feedback from the first generation of massive stars. 
The short-lived sources keep ionizing the baryonic material 
in the halos they form in, as well as their surroundings. How- 
ever, as they turn off, material can cool again and repopulate 
the dark matter hal os. So while the baryon fraction (Fig. 4 in 
IWise & Abelll2008l) fluctuates and decreases at times to as lit- 
tle as 10%, star formation can continue as long as no sustained 
external UV flux sterilizes the halo. The latter case severely 
limits star formation and has been discussed many time in the 
literature (e.g., [Ba bul & Rees 1992; Thoul & Weinberg 1996; 
iKepner et al.|[T999l; IDiikstra et al.ll2004l) . It seems clear then 
from the limited guidance we have from numerical simula- 
tions that most Milky Way satellite halo progenitors experi- 
enced most of their star formation before they are permanently 
ionized. 

Once we have identified satellite galaxies in the simula- 
tion, we must assign magnitudes to them in order to make 
direct comparisons with observations and to account for ob- 
servational completeness effects. This is done using two 
methods. First, we use a halo abundance matching method 
(iKravtsov et al.ll2004al; iBlanton et alJl2008h . Here, luminosi- 
ties are assigned to halos by assuming a one-to-one corre- 
spondence between n{< My), the observed number density 
of galaxies brighter than My, with «(> I'max), the number 
density of simulated halos with maximum circular veloci- 
ties larger than I'max- For the distributi on of magnit udes, we 
use the double-Schechter fit of Blanto net al.l (l2005h for low 
luminosity SDSS galaxies in the g- and r-bands down to 
Mr = -12.375. The Vmax values are taken from the halo catalog 
of a 160 Mpc/h simulation complete down to Vmax ~ 90km/s. 
In order to extrapolate this to lower circular velocities, we 
calculate a power-law fit to the low end of the dyi/dv^nax func- 
tion. The resulting correspondence is shown in Figure [T] for 
the r-, g-, and V-bands (red, green, and black curves). The V 
band magnitudes are calcu lated using the tran sformation V = 
g-0.55(g-r)-0.03 from lSmifli eldl (l2002h . This method 
implicitly ass umes that all galaxie s have average color. Since 
the data from IBlanton et alT(l2005h is not deep enough to map 
onto the dwarf galaxy distribution, we use a power law to ex- 
trapolate the Mv(vniax) relation to lower magnitudes. For the 




Fig. 1 . — The relationship between magnitude and v„,ax for the r-, g—. and 
V— bands using abundance matching (solid red, green and black lines). The 
dashed lines show power law fits to the low-luminosity end. 



V-band, we get 



Mv-51og(/i)= 18.2-2.51og 



Ikm/s/ 



(1) 



When selecting the appropri ate Vmax for assigning a luminos- 
ity, we follow the method of IConrov et alj ( 120061) and choose 
the peak Vmax over the trajectory of the subhalo for subhalos 
that eventually cross the lO^K post-reionization star forming 
threshold. For subhalos that never reach this threshold, we use 
the value of Vmax at Zreion- In both cases, this then corresponds 
roughly to the mass the halo had at the redshift they stopped 
rapidly forming stars. 

The appeal of this method is that we are able to ignore 
much of the poorly understood (and poorly simulated) physics 
of galaxy formation using a statistical method that has been 
shown to, on average, reproduce a wid e variety of observ- 
able properties for more ma ssive galaxies (IConroy et al.l2006l : 
IConroy & Wechsle^ 120091) . as well as some properties of 
dwarf galaxies down to Vmax ~ 50km/ s (|Blanton et al. 2008). 
It is still unclear how this method will fare at lower masses; 
it must break down for small halos once they no longer host 
one galaxy on average. If this transition is sharp, however, 
it may be a reasonable approximation for most of the mass 
range where halos host galaxies. 

As a second approach for assigning magnitudes, we use a 
toy model to predict the star formation rate and stellar mass 
of a satellite combined with the stellar population synthesis 
(SPS) code of lBruzual & Chariot! (I2003h ^ Here, we again as- 
sume that star formation begins when the satellite first crosses 
the mass threshold, Mt, and ends at the reionization time, 
Zreion- During this period, the star formation rate is set by the 
dark matter mass of the subhalo. 



SFR = I ^ [f^'idg^^ 1 Mp 




if Mom >M,, z>z, e 
Otherwise 



(2) 



where fcoidgas is the fraction of cold gas in the halo, and 
a and e are free param eters. This is similar to model IB 
of iKoposov et al.l (120091) . with a couple of key differences. 
First, we impose a hard truncation of star formation at the 
epoch of reionization, something they only consider using 

^ http://www.cida.ve/ bruzual/bc2003 
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their model where stellar mass is a constant fraction of dark 
matter mass. Second, they treat stellar mass as being propor- 
tional Mg^ at some epoch, while we take the total stellar mass 
to be proportional to t he integral of this quantity, as motivated 
by observations (e.g., 'Juneau et al. 2005 ; No eske et al.ll200'7l : 
IZheng et alJl2007tlConroy & Wechs ler 2009). For our model, 
we hold fcoidgas constant and set it equal to the universal 
baryon fraction during the period of active star formation. Our 
two free parameters e and a, set the efficiency and scaling of 
the star formation. We keep e as a free parameter that is used 
to match the luminosity function, but fix a = 2, extrapolated 
from hig her mass galaxies at both low and high redshifts (Fig - 
ure 8 of IConrov & Wechsleij l2009t iDrorv & AlvarezI 120081) . 
The implicit assumption of this model is that subhalos con- 
tain a gas fraction equal to the universal baryon fraction, and 
that this gas exists in one of two phases. At early times the gas 
in low mass halos has a temperature less than T = 8 x 10"'K 
and therefore cannot cool by atomic hydrogen cooling. How- 
ever, once a halo reaches Mt, cooling becomes effective and 
all gas rapidly enters a cold phase where it is able to form 
stars. Reionization, however, rapidly heats all gas in the sub- 
halo to a few times lO'^K, quenching star formation. The more 
massive halos eventually shock heat to such a virial temper- 
ature, which allows the ionized gas to cool again and resume 
star formation. Only 7 subhalos in the via Lactea II simu- 
lation ever reach such a mass. Regardless of mass, all star 
formation is ended once a satellite accretes onto the larger 
host halo, as we assume that this process causes all gas to be 
stripped from the subhalo. Most of the z = subhalos in via 
Lactea //were accreted in the range z = 0-6. This model does 
ignore a large number of physical processes, such as recombi- 
nations and feedback, but much of this can likely be accounted 
for by appropriately setting the constants e and a. In Section 
§5 we discuss how this model compares with previous stud- 
ies. Satellite magnitudes for this model are determined us- 
ing the stellar population synthesis code of lBruzual & CharTol 
(l2003h . For this model, we treat the star forma tion as a series 
of bursts, all with the IMF of Chabrier (l2003h and metallic- 
ity Z = O.OOO4Z0, broadly consistent with the expected level 
of enrichment fro m the earliest, most massive stars tha t pre- 
enrich the halo gas jAbel et al.l2002l:IWise & Abell2008h . Us- 
ing these magnitudes, we tune the e parameter of equation |2] 
so that our model most accurately reproduces the observed lu- 
minosity function of local satellite galaxies (see below). This 
results in typical values for e ~ lOM© /yr. 

As a direct comparison, the resulting average Mvivmux) re- 
lation from these two methods do not differ significantly. 
When considering the populations as a whole, the SPS method 
results in a weaker magnitude trend with v„,ax and slightly re- 
duced magnitude range. On a galaxy by galaxy basis, how- 
ever, this is appreciable scatter, with galaxy magnitudes dif- 
fering by as much as 1.5 dex. 

3. OBSERVATIONAL SAMPLE 

In the next section we present a number of comparisons 
with the observed properties the Milky Way satellites. We 
take as our data set the 11 classical and 12 newly-identified 
SDSS dwarf galaxies within the fiducial mil ky way virial 
radius of 417kpc as presented in Table I of iToUerud et alJ 
(|20(38i). However, we exclude Bootes II from our comparisons 
of the maximum circular velocity, v^ax, because there are only 
poor constraints. When using the SDSS dwarfs, we must ac- 
count for both the sky coverage of SDSS and the detection 
efficiency for each galaxy (.Koposov et al...2008.) . While it is 




Fig. 2. — Distribution of reionization epochs for Milky Way mass halos, 
MpnF = 1.6 — 2 X lO'^/r'M0. The do tted lines and hatche d region indicate 
the 95% limits for the distribution. See lAlvarez et alj 120081) for more details 
on the simulations used to generate this distribution. The dashed red line 
indicates the median value Zmmn = 1 -1 ■ 

Straight-forward to adjust the observational data to account 
for sky coverage, a more subtle issue arises with the depth of 
SDSS. Because the survey is magnitude limited, SDSS is only 
complete in searching for satellite galaxies of absolute mag- 
nitude My out to some radius Rcmp- This completeness depth 
is roughly independent of surfac e brightness and the relation 
is given in lKoposov et al.l (l2008h as 



47r/D 



1/3 



^0(-«M„-fc)/3j^p^^ (3) 



where foRs = 0. 194 is the sky coverage frac tion of SDSS, and 
a = 0.6 and b = 5.23. iTollerud et al.l (120081) also used the via 
Lactea II simulation to calculate a correction due to this ef- 
fect and use it to correct the luminosity function. However, 
this correction is unable to account for the radial or Vmax dis- 
tributions. Consequently, whenever we compare our model 
predictions with observations, we present the sky coverage- 
corrected distribution of observed satellites and a sub-sample 
of the simulated galaxies that is cut by the magnitude Umit of 
equation[3] 

In addition to survey depth, the observational sample 
of Milky Way dwarfs is also limited by a surface bright- 
ness detectabil i ty con straint. This is discussed in detail in 
iKoposov et al.l (120081) who parameterize the detection effi- 
ciency of a satellite according to 



e(M,„Ai) = G 



M,-M,„ 



I im 



(4) 



where G(x) = 0.5erfc(x/V2) and Myji,„, /x//„,, ctm, and de- 
scribe the magnitude and surface brightness limits of observa- 
tions. Using their parameters for fitting current observations, 
we calculate the detection efficiency of our model galaxies us- 
ing the conservative assumption that the light of a satellite is 
uniformly distributed out to the tidal radius of its host dark 
matter subhalo. 

Finally, we must account for the observational uncertain- 
ties in measurements of the magnitudes of the dimmest dwarf 
galaxies. Effectively, the small number of stars used to deter- 
mine the magnitude of a faint dwarf galaxy introduces shot- 
noise effects, creating additional un certainty on the measure- 
ment. This is studied in detail in iMartin et al.l (l2008i) . who 
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report magnitude uncertainties in their Table 1. We include 
these uncertainties when making comparisons between our 
modeled luminosity function and observations. 

4. RESULTS 

4.1. Distribution of reionization epochs 

First, we first investigate the distribution of reionization 
times for MiUcy-Way size halos in Figure |2] The solid line 
of the Figure shows the distribution of reionization epochs 
for the ^500,000 identified halos with mass MfoF = 1 .6 - 2 x 
10'^/!~'Mq from our Gpc simulation. From this plot, we see a 
wide distribution of reionization epochs, peaked around z = 6 
with a tail extending beyond z = 16, fit by an exponential, 
dN/dzreion e"^-*^'"™. Ninety-five percent of the halos are 
reionized in the range z,eion =6.3-11.6 with a median red- 
shift of 7.7. This distr ibution is also quite co nsistent with 
similar predictions from lWeinmann et al.l (120071) . Percolation 
of reionization happens at z = 6.3 and the universe rapidly be- 
comes fully ionized, which is responsible for the sharp cutoff 
in the distribution. This broad distribution indicates that a pre- 
cise constraint on the globally-averaged reionization epoch, 
for example as measured from the optical depth to the cosmic 
microwave background, does not give precise constraints on 
the reionization history of our local galaxy and its progenitors. 
As we will show, determining where the Milky Way itself 
sits in this distribution may be important for understanding its 
satellite population. The reionization epoch of a given halo 
is almost certainly correlated with other properties, including 
the large scale bias and the detailed formation history of the 
halo, and there may be observational clues beyond those pre- 
sented here, to where the Milky Way lies in the distribution. 
We postpone an investigation of these issues to future work. 

4.2. Distribution of satellite populations 

Using the distribution of z,-eion from Figure |2] as a guide, 
we show in Figure |3] the number of subhalos hosting satel- 
lites galaxies, «sats, in the via Lactea II simulation for our 
model. The contours of the lower panel show the full variation 
in Hsats as a function of both Zreion and M^. The dashed lines 
denote constant virial temperatures. The thick lower line rep- 
resents Tvii- = 8 X lO^K, the temperature where HI cooling is 
expected to be effective for gas that has not been photoionized 
(iWise & Abei2008h . The upper fine is lO^K, the mass where 
the halo has sh ock heated to the point w here photoionized gas 
can cool fe.g.. lHaiman & Brvanll2006h . The top panel shows 
Hsats as a function of Zmon when Mi is set such that star for- 
mation begins once the subhalo reaches a virial temperature 
Tyir = 8 X IO-'A'. This line asymptotes to «sats = 7 because we 
allow star formation to resume after reionization for halos that 
shock heat to temperature T^^ = lO^K. Virial temperatures are 
calculated using the relation 



1000 



M= 10" M, 







1800.K 



21 
T+z 



0.3 



1.22 

1/2 



3/2 



■MS) 



0.7 

Ti 



-^]. (5) 



where = 0.57 for T^i, < 1.5 x lO '^K and 1.22 for TVi, > 
1.5 X lO'^K ( iHaiman & Brvanll2006h . The model referred to 
as "all galaxies" in this plot indicates all subhalos capable of 
hosting a galaxy. The galaxies are not defined in terms of a 
magnitude limit, but weather they reached a mass scale (left 
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Fig. 3. — Bottom: The number of surviving satellite galaxies in the via 
Lactea //simulation assuming that subhalos had to grow above a threshold 
mass Mt by Zmon i" order to host a galaxy. The dashed vertical lines show the 
95% confidence intervals for the reionization epoch of Milky Way mass ha- 
los. The dashed red lines denote constant virial temperatures, where 8 X IVr'K 
is the expected temperature for allowing star formation through HI cooling. 
The the light and dark shaded regions show the observational constraints on 
the both the maximal and most like ly distribution in the number of satellite 
galaxies from ITollerud et all J2008I) . The yellow contours denote the 95% 
confidence levels for constraining Zrcion and Mt by matching with the Vmax 
function of Figure [7] Top: The number of satellite galaxies as a function 
of ireion assuming that M, follows the r = 8 X IQ^K curve. The dotted lines 
show the 95% distribution for Zmon for Milky Way mass halos, and the yellow 
lines show the 95% confidence range in Zrcion at this Mt. Again, the light and 
dark cyan regions show full and most likely range for the number of Milky 
Way satellite galaxies. 

as a variable in this model) above which star formation was al- 
lowed to proceed. As a consequence, many such galaxies are 
expected to be extremely low luminosity, well below current 
detection limits. Additionally, the model is simply a function 
of the subhalo mass at times prior to accretion and does not 
track effects such as tidal disruption. 

The most striking feature of this plot is the size of the vari- 
ation in Hsats, which can differ by roughly an order of mag- 
nitude at fixed T(Mt) as Zreion is varied within the 95% distri- 
bution or as Ml changes from Tyir =10^- lO'^K. For a con- 
stant Ml, the number of satellites has a roughly exponential 
dependence on Zi-eion, «sats oc e"" ''5&cio„ Again, there is an im- 
posed minimum value of nsats = 7 for any value of Zreion, set 
by the assumption that photoionized gas can cool and form 
stars when a halo shock heats to Tvi, = lO^K; seven of the 
via Lactea II subhalos reach this temperature at some point in 
their histories. 

The shaded region of Figure[3]shows the observational con- 
str aints on the number of satellite galaxies from the work 
of ITollerud et"an (l2008h . The lighter regions shows the ex- 
treme case of this analysis, where the Milky Way has between 
23 (currently observed) and 2261 satellites galaxies, while 
the darker region shows what their work considers the more 
likely prediction of 300-600 satellites. When all this data is 
combined, the consistency is rather surprising. The region 
spanned by the upper and lower 95th percentile for Zieion of 
Milky Way mass halos, as well as the Mi w IO^Mq/i"' from 
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lAbel et al.l ( 120021) and the Tyii- = 8 x lO ^K falls entirely within 
the constraints of lTollerud et alj ( 120081) . It is important to note 
that these results rely on simulations using completely differ- 
ent sets of physics. The limits on Mi are set from hydrody- 
namical simulations as discussed in ^2.21 while the Zmon and 
«sats limits come from N-body simulations that model only 
collisionless physics. 

One caveat that must be kept in mind when interpreting Fig- 
ure [3] is that we have assumed the particular subhalo popula- 
tion of via Lactea II is representative of a typical Milky Way 
mass halo. While only a handful of such ultra-high resolu- 
tion simulations have been conducted, it is akeady apparent 
that there is a wide distribution in the number of subhalos in 
halos of similar mas s. In particular, currently the three most 
well resolved halos JPiemand et al.ir2008l: IStadel et al.ll2008l: 
ISpringel et al.ll20()8l) contain a factor of 1 .5-2 more subhalos 
than via Lactea II at a fixed mass threshold, and it is estimated 
that via Lactea II is among the 15% of objects with similar 
mass that have so few subhalos (llshivama et al.ll2008b . If this 
is the case, the number of satellites predicted in Figure |3] is 
potentially a factor of two too low for a typical Milky Way 
mass halo. However, it is unknown exactly where in the rela- 
tively wide distribution the Milky Way lies, particularly since 
the number of subhalos has been shown to corr elate strongly 
with halo concentration and formation history dZentner et al.l 
l2005h . In the remainder of this work, we assume that the sub- 
halos in via Lactea II are representative of the Milky Way, but 
this distribution, and the possible bias, should be kept in mind 
when detailed numerical results are given. 

4.3. Luminosity Function 

While Figure |3] shows that the total number of subhalos 
hosting satellite galaxies may be strongly dependent on the 
time of reionization, it is necessary to understand the proper- 
ties of these affected halos, i.e., are they all low mass objects 
that we expect to host low-luminosity galaxies, or do they fill 
a larger range in satellite parameter space? In order to quan- 
tify the expected impact on observations, we must first impose 
the relevant observational cuts on our satellite distribution. 
For each subhalo, we calculate rs„„, the distance from a point 
8kpc from the center of via Lactea II. Figure|4]shows this dis- 
tribution as a function of magnitude for the model Zieion = 8, 
Mt = 3 X 10^/z~'Mq. The open red circles show magnitudes 
calculated using the abundance matching method (equation 
[T]i, and the filled green triangles use the population synthesis 
model. We then impose the cut defined by equation |3] above, 
shown as the black line. Because we expect this subset to 
best match the observational sample, this cut is imposed for 
all subsequent comparisons. While only affecting about 20% 
of our satellites, objects as bright as My = -7 are cut. The 
distributions of the Milky Way dwarfs are overplotted for ref- 
erence. We also calculate the detection efficiency for each 
satellite galaxy based on its surface brightness according to 
equation|4] All galaxies passing out magnitude cut have e ~ 1 , 
so we do not make any additional cuts based on estimated sur- 
face brightness. 

Again, because the magnitudes set by the abundance match- 
ing method are not directly impacted by Zieion and Mt, the dis- 
tribution of objects in My - r j,,,, space in not strongly impacted 
as these parameters are varied. In particular, adjusting these 
parameters only results in the presence or absence of objects 
with low My as low mass subhalos gain or lose the ability to 
host satellite galaxies. Individual objects will, however, have 
a significant dependence on magnitude in the SPS model be- 
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Fig. 4. — The distributions of satellite galaxy magnitudes as a function of 
distance from the sun in a model with <:reion — 

9andM, = 3x IO'/j-'Mq. The 
open red circles show magnitudes assigned using the abundance matching 
method (equation [T) and the filled green triangles have magnitudes assigned 
using the Bruzual and Chariot stellar population synthesis (SPS) code. The 
cyan stars show the distribution of the observed Milky Way satellites. The 
solid line shows the completeness depth of the SDSS survey as given by 
equation [5] 

cause adjusting these parameters impacts how long star for- 
mation is allowed to proceed for, impacting the amount of 
mass that can be converted into stars. In addition to form- 
ing new satellite galaxies, pushing Zieion to later epochs also 
causes the existing satellites to brighten. 

Figure[5]compares the luminosity functions from our model 
with observations. The thicker lines show magnitudes set by 
the abundance matching method, and the thinner lines by the 
SPS model. For this plot, we have fixed Mt to be set by the 
rvii-(Mt) = 8 X lO^'K relation and varied z,eion = 5, 8, and 12 (red 
dotted, green dashed, and blue dot-dashed lines). The black 
long-dashed line with points shows the observed luminosity 
function, while the cyan region shows the Poisson errors about 
this distribution. 

Overall, both the abundance matching and SPS methods re- 
produce the observed luminosity reasonably well. The level 
of agreement for the abundance matching method in particu- 
lar is rather impressive since the method 1) ignores all baryon 
physics such as star formation, 2) assumes all galaxies have 
average color, and 3) extrapolates number densities down to 
regimes where the method has not been tested and where 
small-scale processes may produce a significant amount of 
scatter in the v^-j^^-My relation. Because the primary effect 
of changing Zreion is to change the total number of satellite 
galaxies, this serves to change the overall normalization of 
the luminosity function while retaining the slope. Although 
not shown here, changing Mt for a fixed Zieion has a similar 
effect. 

The high level of agreement between our SPS model 
and observations comes from tuning the efficiency param- 
eter, e, to match observations independently for all values 
of Zieion- Thus, the three curves in Figure |5] have values 
6 = 0.1,10, and 3OOM0 /yr, which were selected for no phys- 
ical reason other than to match the luminosity function. It 
is, however, interesting to note that the faint-end slope of this 
model almost perfectly matches the observations. The under- 
prediction of luminosities at the bright end can potentially 
be explained by residual, ongoing star formation. At least 
one of the classical Milky Way Dwarfs, Leo I, shows signifi- 
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Fig. 5. — Luminosity functions for observations and model predictions. 
The long-daslied line shows the observed Milky Way satellite luminosity 
function corrected for sky coverage and depth effects, while the cyan swath 
represents the statistical error The red dotted, green dashed, and blue dot- 
dashed lines represent reionization models of varying Zreion — 5,8, and 12, 
respectively. Mt is set using the virial temperature, rvii(Mt) = 8 X lO^K. The 
thicker set of lines shows predicted luminosity functions using an extrapo- 
lated abundance matching method to assign luminosities to the galaxies. The 
thinner set of lines use a stellar population synthesis model to predict the 
luminosities. 



cant evidence - while Fornax shows slight evid ence - for re- 
cent star formatio n after the epoch of accretion (iMateoll 19981 : 
iMateo et"ani2008l) . If we were to allow some such process in 
our 4 brightest objects the model will likely fit the observa- 
tions significantly better 

Finally, we consider the ration of mass to luminosity in the 
top panel of Figure |6] Here, we plot the mass to light ra- 
tio of the via Lactea II subhalos using their virial masses at 
time of accretion and luminosities assigned from the abun- 
dance matching method for z,eion = 8 as red circles. Solid 
circles represent objects within the SDSS magnitude limit, 
rsun < RcompiMv), wherc Rcomp is given by equation[3] Open 
circles are satelUtes outside this limit. This model naturally 
reproduces a wide range of M/L rations, spanning from 10^ 
to 10*", with a clear trend with luminosity, M/L cx L""^^. 
Green triangles represent luminosities calculated using our 
SPS method. This method reproduces a similar trend with 
luminosity but with a larger dispersion. We have included 
measurements from the MiUcy Way as the cyan stars. In or- 
der to model the masses of these objects, we took the Mp 3 
values for the mass within 300kpc published in lStrigari et alj 
(12008^) and converted those to subhalo masses using their pub- 
lished relation M0.3 = lO^M0(M„>/lO''Mo)°■^^ where M,„> is 
the virial mass of the subhalo at the time of accretion. We 
must caution that this relation is expected to be dependent on 
cosmology and ignores all scatter. The observations are re- 
markably well matched by our model with excellent agree- 
ment for all but the most most luminous galaxies. In the 
lower panel of Figure |6] we consider this data in a different 
way by plotting M0.3 as a function of luminosi ty. Here, we 
take the data directly from IStrigari et al.l (l2008h and convert 
the via Lactea II subhalo masses to M0.3 using the above for- 
mula. While the numbers are in general agreement, the abun- 
dance matching model (red circles) shows a clear trend with 
luminosity, M0.3 cxL° '\ as opposed to the observations (cyan 
stars) which indicate more of a common mass scale. The trend 
in our model results directly from the abundance matching 
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Fig. 6. — Top Panel: The ratio of subhalo mass at the time of accretion 
to luminosity for our abundance matching model (red circles) with z,cion = 8. 
Filled circle represent halos that are within the SDSS completeness radius, 
equation[3] while open circles are outside this distance. Cyan stars represent 
Milky Way satellites. In order to make this comparison, we used the values 
for ma ss within 300 kpc at the present epoch. Mo 3, published in Strigari et al] 
1200 8ft and converted them to subhalo masses using their published relation 
for the average Mdm(Mo 3) calibrated to N-body simulations. Bottom Panel: 
The relation between Mo 3 and luminosity for our abundance matching model 
and observations. Here, we have converted the mas ses of the via Lactea 
II subhalo to Mq 3 values using the formula provided byl Stri g ari et alJ (20081) . 

method, which assigns luminosities to subhalos satellite based 
on Vtnax at the time of accretion. The addition of scatter into 
either the Lv{vm-dx) or M().3(M,,,>) relations, which we expect 
at these low mass scales, can help to flatten this trend slightly 
and bring it more in line with observations. The SPS model 
(green triangles) also produces a similar trend with luminosity 
albeit with a significantly larger scatter, making the slope of 
the Mo j,(Lv) relation consistent with zero. Prev ious studies of 
simulations have observed a similar trend (e.g.. lMaccio' et af] 
'20Q8U Li et alj20 08VKopo sov et alj2009t) using semi-analytic 
modeling of galaxies and/or subhalo distributions. 

4.4. Circular Velocity and Radial Distributions 

We next consider the mass distribution of the satellite galax- 
ies hosting halos. Figure |7] shows the changes in the Vmax dis- 
tribution for satellites as Zi-eion is varied, given a threshold mass 
of rvir(Mt) = 8 X lO^K as in Figure|5] Here, the solid black Une 
shows the distribution for all via Lactea II subhalos, while the 
red dotted, green dashed, and blue dot-dashed show the distri- 
butions from our model for three values of Zreion = 5, 8, and 12. 
Because the abundance matching method was more success- 
ful than the SPS model in reproducing the observed luminos- 
ity function without the need to tune any parameters, we only 
include satellites that pass the radial cut of equation [3] using 
the abundance matching criteria. As can be seen, an earlier 
Zreion supprcsscs the distribution of subhalos with all values of 
Vmax, although the effect is more pronounced for low mass ha- 
los. Still, this suppression is present even for v^ax ^ 20km/s, 
where most of the classical dwarfs live. This indicates that 
Zreion Can effect satellite galaxies of all masses and luminosi- 
ties. 
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Fig. 7. — The peak circular velocity functions for subhalos hosting satel- 
lite galaxies. The solid line shows the velocity function for all via Lactea 
II subhalos. The red dotted, green dashed, and blue dash-dotted lines show 
the distribution for subhalos hosting satellite galaxies for i,cion =5,8, and 1 2 
with rvir(M) = 3 X 10'*K, as in Figure [5] The long dashed black line shows 
observed distribution for Milky Way satellites, corrected for sky coverage 
and detection efficiency. The cyan bands show combined Monte Carlo and 
statistical errors. 



The long-dashed black hne with data points again rep- 
resents the observations. The v^^x values for the satel- 
lites, including error s, were calculated using the method of 
iStrigari et a l. (2007b a) using kinematic data taken from the 
hterature (jWalker et al. 2007,' Simon & Geha 2007). The line 
was calculated using the 22 observed satellites and correct- 
ing them for SPSS sky coverage and detection efficiency 
jKoposov et al.ll2008h . The cyan region denotes errors on 
this curve and were calculated using a Monte Carlo approach. 
In this approach, published errors are used where possible; 
where no robust errors are published, the average error distri- 
bution is mapped onto the remaining SDSS dwarfs. While this 
process does not produce uncorrected error bars, it should be 
significantly more robust than simply using statistical uncer- 
tainties. We should also note that, because the reconstruction 
of Vmax from observations gives a very strong lower bound but 
only a weak upper bound, the errors on the lowest points are 
probably underestimated because very few satellites will scat- 
ter into this bin as we create a Monte Carlo representation of 
the distribution. These systematic errors were combined with 
statistical errors assuming a Poisson distribution. 

Figure[8]further explores the impact of varying Zieion on the 
properties of the satellites and the subhalo hosts by consider- 
ing the radial distribution within the halo. The lines represent 
the same populations as in Figure [7] Here, the errors on the 
observations are purely Poisson. Again, there is strong trend 
for an early reionization epoch to suppress the abundance of 
satellite galaxies at all radii, but, in part because the data is 
rather noisy, it remains easy to match the model to the obser- 
vations. 

Finally, using our magnitude-limited sample of satellites 
from the abundance matching method, we can attempt to put 
constraints on Zieion and M^. We calculate the errors be- 
tween the Vmax distribution for observations and our models 
and plot the 95% confidence levels as the yellow contour in 
Figure |3] We use the Vmax distribution to generate this con- 
straint because it is less sensitive to the modeling of the mag- 
nitudes than the luminosity function is, with the modeling 
only coming into play when making the SDSS magnitude- 



FlG. 8. — The radial distribution subhalos hosting satellite galaxies in via 
Lactea II. The lines represent the same populations as in Figure|7] 

limit cut. Because the Monte Carlo method underestimates 
the error on the lowest v,nax point in Figure |7] we exclude it 
from our constraints. The general insensitivity of our con- 
straints to the value of Mt comes from our use of VmnxiZraon) 
in setting the magnitude of a galaxy. While lowering the 
threshold for star formation from ~ 10*^ to ~ \(fh~^M.Q sig- 
nificantly increases the number of subhalos hosting galaxies, 
most of these objects are quite dim, with My - 4, and are 
therefore cut from our observable sample (see Figure |4|i. If 
we impose the condition that galaxy formation begins when 
a halo reaches virial temperature T^i^iMx) = 3 x lO^K we can 
constrain the time of reionization to Zreion = 8^,, with 540^34q 
satellites. Note that these errors assume that the dark matter 
substructure of the Milky Way is identical to that of the via 
Lactea II halo; if the Milky Way is more typical for its mass 
it may have a larger total number of observable satellites. A 
further success of this model is that, when we constrain Zmon 
using the v^ax function, we naturally reproduce both the ob- 
served luminosity function and radial distribution, as shown 
by the green dashed lines in Figures |5] and [8] We should, 
however, caution again that much of this depends on the par- 
ticular realization that is the via Lactea II halo, and that more 
statistics will be necessary for a more robust prediction. 

5. COMPARISON WITH PREVIOUS WORK 

Our result that the number of satellite galaxies is strongly 
dependent on the redshift of reionization i s at odds wi t h som e 
previous studies, inclu ding t he work of ISomervilld ( 120021) . 
iKravtsov et all ( l2004bl) . and lOrban et alj ( l2008h . although 
these studies primarily focused on the classical satellite galax- 
ies in the Milky Way. The primary difference in our models 
rests on the assumption of how reionization effects the pres- 
ence of cold gas. These previous studies have used the model 
of Gnedin (2000) to calculate the amount of cold gas in a halo 
of mass M. They calculate a filtering mass, Mf{z), the mass of 
a halo that looses half its baryons compared to the universal 
baryon fraction in the presence of a photoionizing UV back- 
ground. This mass is related to the baryon fraction of a halo 
of mass M via the relation 



/gas(M,z): 



/baryo 



[I+0.26M/(z)/M]3 



(6) 



where the details for calculating the filtering mass, Mf, are 
given in [Gnedia (200Q) and Appendix B of IKravtsov et alj 
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(I2004bl) . 

The crucial point in using the values of /gas from iGnedinI 
( I2OOOI) lies in how these fractions are related to the amount of 
cold gas available for star formation. There are two bracket- 
ing possibilities: either the gas is spread out over the entire 
halo and is hot, in which case the star formation rate is zero, 
or, there is a small clump of cold gas in the center, with a mass 
given by /gas^haio, that can form stars. Our interpretation, in 
which the gas is assumed to be hot for /gas ^ 1, is consistent 
with the first possibility, while these previous studies have im- 
plicitly assumed that all the gas is cold. 

There are seve ral reasons to favor the hot gas scenario. 
Most importantlv. iGnedinI ( 120001) did not distinguish between 
hot and cold gas when calculating the gas fraction. The as- 
sumption that the gas is cold is only valid for halos that had 
collapsed before the reionization epoch, and were then sub- 
ject to photo evaporation due to the UV background. Such 
cold gas would only survive for one photo evapor ation time, 
which , according to the numerical simulations of llliev et alj 
(l2005h . is likely to be less than 500 Myr, corresponding to 
Az < 2 at z < 6. After such photo evaporation, the halo 
could only accrete hot gas, which would be unlikely to form 
stars given its low density and long cooling time. Moreover, 
due to the exponential growth of the abundance of halos with 
masses ~ 10** M© during reionization, it is likely that most 
halos present after reionization were only just collapsing, and 
thus did not have any cold gas capable of forming stars in the 
first place. We therefore interpret the gas fractions /gas <C 1 
reported by lGnedinl (12000!) after reionization as corresponding 
to hot gas when /gas <C Vli,/Vl„,. 

Another related caveat in using /gas to model the depen- 
dence of star formation history on reionization is that the 
value of /gas reported by IGnedinI (|2000) is a mean value, for 
all halos in the box, regardless of whether they are in ionized 
regions or not. In reality, those halos that were reionized ear- 
lier have lower gas mass fractions than those that were reion- 
ized later. Using the average value of /gas for all halos under- 
estimates the sharpness of the transition for halos that were 
ionized at a given time. Because our reionization redshift is 
defined for a given halo (as opposed to a universal time for the 
universe), we expect the transition in gas fraction to be much 
sharper than that given by averaging ov er all halos, eac h with 
differing values of Zi-eion, as was done bv lGnedinI (120001) . 

The difference between these scenarios — where fgas rep- 
resents the fraction of hot in a halo consistent with our abun- 
dance matching and SPS models as opposed to the fraction 
of cold gas available for star formation — is illustrated by 
Figure |9l which shows the evolution of faasiM,z), the ma- 
terial available for star formation, with redshift. The upper 
and lower pairs of lines in the figure represent halos with 
mass M = lO'^/^'M© and 1O''/!"'M0. The soHd and dashed 
lines represent models where the universe reionizes at differ- 
ent epochs, z = 5 and 12, respectively. At high redshift, ha- 
los have a gas fraction equal to the universal baryon fraction. 
As time evolves, the /gas of the average halo decreases rather 
slowly, due to photoionization heating as halos begin to be ex- 
posed to the UV background. Changing the redshift of reion- 
ization affects the rate of transition from the cold to the hot 
phase, but the overall shape of the transition is preserved. As 
we emphasized above, because the average over all halos was 
used to calculate the filtering mass, this is likely an accurate 
description of how reionization effects the gas content of an 
average halo of mass M given a universal redshift epoch. We 



lU" 




10"- 








- 


10-" 










3 10"' 




















10-'° 










10-'^ 
10-" 


_ M = 
. M = 


10*Ms 
lO'Ms 







15 10 5 



z 

Fig. 9. — The fraction of mass in the form of cold gas as a function of 
time. The solid and dotted Hnes are calculated using universal reioniza- 
tion epochs z = 5 and 12, respectively. T he lowe r (blac k) and upper (red) 
pairs of curves show the predictions from IGnedinI I2000f) for halos of mass 
Mom = 10*/i"'Mq and lO^/r'M©, respectively. The vertical lines represent 
the model used in this work, where the reionization time of the subhalo is 
Zicion = 6.5, 11.5 for the solid and dashed lines. The horizontal dotted line 
shows the universal baryon fraction. 

contrast this with the model used for this work, shown by the 
vertical lines indicating a sharp cutoff in the presence of cold 
gas at the time at which the halo reionizes is mass indepen- 
dent. While we acknowledge that we are ignoring effects such 
as atomic recombination and additional cooling for high mass 
objects that should not quite reduce the cold gas fraction to 
zero after reionization, we expect that such a sharp transition 
from the cold to hot phase more accurately describes the evo- 
lution of an individual halo. 

The difference in the resulting star formation histories in the 
two scenarios discussed above is dramatic. To show this ex- 
plicitly, we have re-run our population synthesis model using 
the more gradual star-formation squelching model of equa- 
tion |6] instead of an abrupt squelching for setting fcoUgas in 
equation |2l although we still assume an abrupt end to star 
formation when the subhalo is accreted onto the host halo. 
The result, shown in Figure [Tol is that the luminosity func- 
tion of the satellite galaxies becomes largely independent of 
the reionization epoch, in agreement with these previous stud- 
ies. Here, the thin lines represent our model with an instanta- 
neous gas heating, while the thick lines show the gradual turn 
off from equation |6l The red dotted, green dashed, and blue 
dot-dashed show the effect of changing the reionization time, 
Zreion = 5,8,12. We should note that, unlike in Figure|5] we did 
not tune the parameters a and e of equation[2]to reproduce the 
observational sample at all values of z,-eion, but just set the pa- 
rameters to fit the Zieion =9.6 model. Additionally, this plot 
does not include any observational magnitude cuts, equation 
[3] While the instantaneous-squelching model shows dramat- 
ically different luminosity functions as z,eion is changed, the 
gradual cutoff model is remarkably stable. This is because, 
even though the rate at which gas heats varies with reioniza- 
tion time, the gradual turnoff causes roughly the same amount 
of stellar mass to form in an average halo regardless of the 
reionization time. The average stellar mass, {\og{Msteiiar)), 
changes by less than 5% as Zreion varies from 12 to 5. For our 
rapidly truncating model, however, the amount of gas con- 
verted into stars clearly depends on the time of reionization. 
This raises an important point, in that while there is some dis- 
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Fig . 10. — The luminosity functions from our SPS model for both the case 
of an instantaneous gas heating (thin lines) and a more gradual turnoff as 
given by equation [6] (thick lines). Colors and line-styles represent different 
epochs of reionization as in Figure|5] 

agreement as to the impact the time of reionization has on the 
satellite galaxy population of a halo, the rate at which star for- 
mation is squelched can also have a significant impact on the 
population. 

In order to understand why the luminosity function does 
not vary with Zreion in this model, it is necessary to look at 
the magnitudes of individual galaxies. The most massive sub- 
halos in via Lactea II are larger than the filtering mass, Mf, 
and are therefore unaffected by changes in the reionization 
epoch. In order to understand why the dim end of the lu- 
minosity function has such a weak dependenc e on the reion- 
ization epoch, we note that, as pointed out in iKoposov et alj 
(|200^, halos in this model with a gradual turnoff of star for- 
mation create a large number of stars after Zieion- Additionally, 
as seen in Figure|9] there is a transition epoch around z = 3 be- 
tween suppression and enhancement of /gas for early and late 
reionization times. At early epochs, an early Zreion suppresses 
the amount of gas in a halo of a given mass relative to a late 
reionization because of the additional energy input to the sys- 
tem. However, at late times, the early Zmon actually causes an 
enhancement relative to later reionization because the expan- 
sion of the universe since reionization causes adiabatic cool- 
ing. This transition is roughly independent of mass. Thus, 
there are two regimes; Subhalos that accrete onto the main 
halo early, before z = 3, must get dimmer as Zj-eion increases, 
because star formation is suppressed at all epochs. However, 
objects that accrete recently pass through a phase where ear- 
lier reionization enhances star formation, potentially allowing 
such satellites to brighten with. 

This trend between star formation and time is seen in Figure 
nn which plots A(Mv/), the difference in magnitude for indi- 
vidual satellite galaxies between early (Zreion =12) and late 
(Zreion = 5) rcionization, as a function of accretion time within 
the context of a gradual star formation turnoff. Positive val- 
ues represent satellites that brighten with early reionization, 
while negative values represent objects that get dimmer with 
early reionization. There is a clear trend with ZaccreiUm'- Prior to 
z = 3, all objects have negative values for A(Mi/), as expected 
from Figure|9] while objects that accreted more recently have 
both positive an negative values, with the most positive val- 
ues (strongest brightening from earlier reionization) coming 
from the most recently accreted satellites. Thus, while Zreion 
does effect the magnitude of a given satellite in the context of 




Fig. 11. — The change in magnitude for satellite galaxies between late 
(z = 12) and early (z = 5) Zmon using the model motivated by|Gnedin 1 200(|), 
where the amount of cold gas is given by equation|6] Galaxies that have been 
accreted more recently are brighter for earlier i,cion due to the boost in star 
formation from adiabatic cooling. The opposite trend is present for galaxies 
that are accreted earlier. 

a model with a gradual star formation turnoff, there are two 
competing effects that cause the luminosity function to be un- 
changed. T his explains the di fferences between our results 
and those of lSomervilld ( 120021) who concluded that Zieion has 
no effect on the satellite galaxy luminosity function. Hydro- 
dynamic al simulations also tend t o dis favor this scenario. A s 
shown in I Abel & Haehneltl (Il999h and lSokasian etal.1 (l2002h . 
quasar activity around z 3 - 4 would have resulted in an 
epoch of Hell reionization, so mething not considered in the 
calculations of iGnediiil (120001) . This would have resulted in 
additional heating of the gas at these epochs, likely destroy- 
ing both the enhanced recent star formation activity due to 
early reionization and independence of the satellite galaxy lu- 
minosity function on the epoch of reionization. 

We also reach differe nt conclusions than the recent work 
of IKoposov et al.l (l2009h as to whether an abrupt or a gradual 
star formation truncation more accurately matches observa- 
tions. In particular, they note a bimodality in the luminosity 
function for an abrupt cutoff to star formation. Their model 
differs from ours, however, in the key respect that they as- 
sume a halo converts a constant fraction of its gas mass into 
stars, resulting in stellar mass cx Mdm, where Mqm is the 
dark matter mass at either reionization or z = 0, depending on 
the mass of the halo. This results in a low-luminosity popula- 
tion of galaxies unable to form stars after reionization, and a 
higher luminosity population that can form stars. If we adopt 
such a model and apply it to the via Lactea II subhalos we 
reproduce a similar set of populations, with a population of 
dim satellites in a very narrow luminosity range separated 
from a brighter population made up of the 7 most massive 
subhalos that are able to sustai n ongoing star f ormati on af- 
ter the reionization epoch. As IKoposov et al.l (120091) note, 
this is a poor match to the observations. However, in this 
work we adopted a more physically motivated model where 
the star formation rate depends on the halo mass, resulting in 
(X i M^udt (e.g., Juneau et al. 2005; Noeske et al. 20021 
IZheng et al. 2007; Conroy & Wechsler 2009), where we inte- 
grate over the time between when the halo first crosses Mt and 
the reionization epoch. Changing the parameter a directly 
allows us to adjust the slope of the luminosity function and 
merge the two populations into a single, continuous distribu- 
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tion. In particular, the value a = 2, measured from high mass 
objects jConrov & Wechsleij 120091: iDrorv & Alvarej 120081) 
provides a close fit to the data. Figure |5] Thus, we conclude 
that a sharp cutoff of star formation can closely match the ob- 
servational data. 

6. DISCUSSION AND CONCLUSIONS 
6.1. Additional Constraints and Physical Processes 

Here we briefly discuss additional physics processes not ad- 
dressed earlier. One observable constraint on our model is the 
presence of fossil galaxies - objects that have had no star for- 
mation since reionization. Current observational constraints 
indicate that all satellites brighter than My ^ -8 have had 
some sort of star for mation in the past 8 Gyr (see the review of 
iTolstov et al.ll2009l and references therein), but that dimmer 
objects such as Herc ules may have had all their stars formed 
during a single event jKoch et al.ll2008h . Our base model pre- 
dicts a large number of fossil galaxies. As noted in section 
2.2, there are only seven "non-fossil" galaxies, whose mass 
reaches our post-reionization threshold of 10''K. The number 
of fossil galaxies in our model, however, is highly sensitive to 
the exact value of this threshold, while the overall abundance 
of satellite galaxies is not. The full impact will be explored 
further in a companion paper, but lowering the second cool- 
ing threshold from lO^K to a few lO'^K can match the fossil 
constraints while not significantly impacting the predicted ob- 
served luminosity function. 

Another limitation of our model as presented is our in- 
complete treatment of tidal disruptions on the satellite galaxy 
population. While the high resolution via Lactea II simula- 
tion automatically models the impact of tidal stripping on the 
dark matter halo, we have ignored the further impact on the 
stellar population. This is equivalent to the assumption that, 
once formed, a satellite galaxy remains undisrupted unless 
its host dark matter halos is completely disrupted, merging 
with the central object or dropping bel ow the resolution limit 
of t he simulation. Recent studies of iMaccio' et al.l (l2009h 
and iKoposov et al.l ( l2009i) have indicated that this may have 
a significant impact on the z = satellite l uminosity func- 
tion. T o te st this, we consider t he criteria of iTaylor & Babull 
(120041) and lMaccio' et al.l (l2009h for disruption, assuming that 
a satellite galaxy is destroyed when the tidal mass of its host 
dark matter subhalo is less than Mbind = M{< fdisfbmd), where 
rhi„d is the radius of the subhalo, at infall within which the to- 
tal energy is positive and i s a fre e parameter that they set 
to fdis = 0.1. ITaylor & BabuH (12004 calculate rbM = 0.11 r„ 
where is the scale radius for an NFW profile. In order to 
estimate Mbind, we adopt the mean relation for concentration 
given mass an d redshift, c{M z) = 9 (M„v/M.^)~" '^ in a ACDM 
universe from lBullock et al] ( l2001h . Using this, we see that, in 
our best fit model with Zreion = 8, roughly 20% of our satellite 
galaxies are effected at all magnitudes. A further investiga- 
tion of this process is needed, but it appears that there may be 
a significant degenracy between tidal disruption and reioniza- 
tion redshift when attempting to match the luminosity func- 
tion. 

6.2. Future Work 

The analysis in this paper opens the door for a significant 
number of future studies. In particular, better statistics from 
both simulations and observations are necessary to distinguish 
between the models discussed here and confirm the result that 
Zreion has a significant impact on the satellite galaxy population 



of a Milky Way mass halo. 

On the simulation side, we need further analysis of both 
large scale and high mass-resolution simulations. By pro- 
viding excellent statistics, large scale simulations of regions 
from 100-1,000 /!"'Mpc will yield a significant amount of 
information about a wide range of properties of Milky Way 
mass halos, including environmental effects on both the reion- 
ization time and subhalo population as well as any correla- 
tions between the two. Clearly, there should be some cor- 
relation between clustering or the proximity to groups and 
clusters with z,-eion since such massive objects are the earli- 
est sources of photoionizing radiation. Addi t ionally, studies 
of halo assembly bias (e.g., iGao et af] 120051 : IWechsler et alJ 
120061: iHahn et al.l2008l) have shown correlations between halo 
assembly history, clustering, and substructure population. It 
is therefore likely that these properties will also correlate with 
the reionization time. Understanding such relations could al- 
low us to measure the environment of the local group to both 
more robustly understand the subhalo population of the Milky 
Way and provide an alternate estimate for its reionization his- 
tory. Such measurements could help confirm or refute the 
model presented here. 

So far we have only applied our model to simulations of 
two single halos - via Lactea II and the original Via Lactea. 
In both simulations, we are able to reproduce the satellite lu- 
minosity and distribution functions to similar accuracy, albeit 
withzreion = 8 for via Lactea II and Zi-eion =10 for the origi- 
nal simulation. By studying additional, more highly resolved 
halos, we can both strengthen our predictions and accurately 
quantify the expected scatter in satellite galaxy population 
given Zreion or an expected «sats for our model. Higher resolu- 
tion will also allow us to accurately measure the distribution 
of mass down to scales of 300 kpc, easing comparisons with 
observations as in Figure |6l Indeed, these scales are just be- 
low the resolution limit in the mos t recent simulations (e.g., 
iDiemand et d]l2008l:[Spnngel et al.l 2008: Stadel et al. 2008). 
Finally, new high resolution simulations of the detailed hy- 
drodynamics of star formation in the early universe can give 
us a better handle on how rapidly we expect reionization to 
truncate star formation by more accurately modeling the rate 
at which the gas in small halos is heated. Such simulations 
could also give independent predictions for the values of a 
and e in equation^ providing another way to test our model. 

There are also a number of additional observational mea- 
surements that would help in understanding the validity of the 
models presented in this work. Most directly, deeper surveys 
will aid in quantifying the abundance of substructure by both 
extending the luminosity function to dimmer galaxies within 
the Milky Way, and providing measurements of other systems 
such as M 31. Upcoming surveys such as PES, PanSTARRS, 
and LSST ( lAbbott et alfcoOSHKaiser etai1l200ailvezic et alj 
l2008l) will map the distribution of galaxies at more than 3 dex 
deeper than SDSS over the entire sky, potentially discover- 
ing hundreds of new satellites. These improved observations 
may also give us a better handle on mass measurements of 
the satellite galaxies. The depth will also allow us to poten- 
tially probe the bright end of the satellite luminosity functions 
for thousands of Milky Way mass galaxies at distances out to 
~ 6QMpc. As with additional simulations, the statistics pro- 
vided by these observations will provide significant discrimi- 
nating power when applied to the models of this work. 

Additionally, there is hope that detailed studies of the star 
formation history of the local satellite galaxies can help by 
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directly measuring both the reionization history and rate at 
which reionization quenches star formation. Because the 
gradual star formation truncation given by equation |6]results 
in a significant amount of stars forming after Zi-eion, detailed 
modeling of star formation may be able to discriminate be- 
tween models with rapid and gradual truncation of the star 
formation rate. Unfortunately for such modeling, any star- 
burst activity occurring after a satellite accretes onto the host 
halo would greatly confuse the results. Because of this, we 
will likely have to concentrate on the dim end of the satellite 
luminosity function, restricting such studies to the Milky Way 
and potentially M3 1 . 

6.3. Summary 

Following the recent work of lAlvarez et al.l ( l2008h . we pre- 
dict a broad range of reionization times for Milky Way mass 
halos, ranging from Zreion ~ 6-12. We find that the time 
of reionization can have a significant impact on the satellite 
galaxy population of a Milky Way halo. We investigate pre- 
dictions for a simple model where, in order to cool gas and 
form stars, a subhalo must reach a threshold mass, Mt, by 
Zi-eion, the time it reionizes. This model predicts a strong de- 
pendence of the satellite galaxy population on Zieion; we find 
that the number of satellites can vary by an order of magnitude 
for a fixed Mt. This result is in contention with a number of 
previous studies which have shown minimal impact of Zi-eion 
on the satellite population. The differing results are due to 
differing assumptions about the rate at which the UV back- 
ground squelches star formation. Previous studies have used 
a gradual tr ansition from th e cold to hot gas phases, based the 
the work of lGnediiil (120001) to predict the amount of cold gas 
available for star formation retained by a halo during the pro- 
cess of reionization. Instead, we interpret this calculation to 
be more indicative of the total amount of gas in the halo, and 
assume that it is rapidly heated to a hot phase so that star for- 
mation is very quickly stopped. Additionally, if the heating 
process causes a slow star formation truncation, we believe 
that the quasar Hell reionization at z ~ 4 - 3 will alter the star 
formation history such that Zi-eion will more strongly impact 
the satellite luminosity function than previous studies have 
shown. While these two interpretations of rapid and grad- 
ual gas heating bracket the most extreme interpretations, their 
discrepancy indicates that further study into the exact heat- 
ing rate is necessary since they predict completely different 
dependencies on Zieion- 

Assuming a rapid heating of the gas from the photoionizing 
background, the reionization redshift Zreion impacts the ability 
of subhalos of nearly all masses and radial distribution to host 
satellite galaxies, with a strong impact on the satellite galaxy 
luminosity function. In spite of the simplicity of the model, 
the results are consistent with both observations and detailed 
hydrodynamical simulations of stellar formation. Extrapolat- 
ing the observed bright end relation between Vmax and lumi- 



nosity down to satellite-sized objects, we are able to closely 
reproduce the luminosity function of the Milky Way tuning 
only the parameters of reionization. This is an additional in- 
dication of the robustness of the abundance matching method 
for assigning luminosities to dark matter halos. 

Because it is minimally impacted by methods for assign- 
ing magnitudes to galaxies, we can use observations of the 
satellite Vmax function to place constraints on the reionization 
epoch of the Milky Way. This observed distribution is best 
recovered for a reionization Zieion = 8^2, slightly more recent 
than the "instan taneous" value of 11. 2 from the WMAP5 data 
jKomatsu et al.l l20Q8). This v alue is also in excellent agree- 
ment with the predictions from lWeinmann et al.l (l2007h for the 
reionization redshift of a MiUcy-Way mass halo. This model 
predicts that the Milky Way should host roughly 540 satellite 
galaxies. Such a value for Zi-eion also produces a good agree- 
ment for the radial distribution of satellites and the luminosity 
function when the abundance matching method is used. It 
must still be cautioned, however, that this result depends on 
assumption that the via Lactea II simulation is representative 
of the dark matter distribution in the Milky Way. We get sim- 
ilar agreement when we adopt a model where SFR cx M^^j^, 
although it is important to note that this model has a tunable 
parameter Still, it is able to reproduce the slope of the dim 
end of the luminosity function almost exactly and some late 
time star formation, consistent with observations, will help 
relieve the (slight) tension at the bright end. Thi s is consis- 
tent with the recent work of Kopos ov et al.l (120091) . Our work 
lends further strength to the growing body of research that 
suggests that there really is no "missing satellite" problem for 
the Milky Way and that the next generation of surveys may 
allow us to understand the entire population of local satellite 
galaxies. 
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